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Abstract 

The fact that several important effects in multiparticle dynamics, on which 
QCD has not yet satisfactory predictions, have been interpreted in terms of 
the validity of negative binomial (Pascal) regularity and related clan proper- 
ties at the level of simpler substructures, raises intriguing questions on clan 
properties in all classes of collisions, the main one being whether clans are 
observable objects or merely a mathematical concept. We approach this prob- 
lem by studying clan masses and rapidity distributions in each substructure 
for e + e~ annihilation and hh collisions, and find that such properties can in- 
deed characterise the different components. These results support the idea 
that clans could be observable, a challenging problem for future experiments. 
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1 Introduction 



As well known, the concept of clan has been introduced in the eighties Jl] in order to 
interpret the occurrence of the approximate NB (Pascal) regularity in final charged parti- 
cles multiplicity distributions (MD's) of the full sample of events both in full phase space 
and in restricted rapidity intervals in all classes of collisions in the GeV region. 

Clans were defined as group of particles of common ancestor with at least one particle 
per clan; by assumption no correlations exist among clans. Accordingly, the production 
process was understood as a two-step process: to independently produced clans in the 
first step (they are Poissonianly distributed), it follows a second step in which each clan 
decays into final particles with a logarithmic distribution. The average number of clans, 
N, in rapidity intervals at fixed cm. energy and at various cm. energies, as well as the 
average number of particles per clan, n c , characterised fully the multiplicity distribution 
properties in various classes of high energy collisions. These two parameters are linked to 
standard NB (Pascal) distribution parameters n (the average charged multiplicity) and k 
(1/k = D 2 /n 2 — 1/n, D being the dispersion of the MD) by two non-trivial relations 



We learned [0], within the just recalled elementary interpretation, that N in e + e~ 
annihilation is larger at the same cm. energy than in hadron-hadron collisions, and that 
n c is much smaller in the former than in the latter case. In deep inelastic scattering 
(DIS), the situation turned out to be intermediate between the previous two: N behaves 
as in hadron-hadron collisions but n c as in e + e~ annihilation. In addition, in all classes of 
collisions, clans are larger (they contain more particles) in central rapidity intervals than 
in peripheral ones. 

These remarkable properties of clans were obtained in a quite simple framework and 
suggested an interesting clan picture at parton level by using generalised local parton- 
hadron duality (GLPHD) H. It was found that partonic clans behave as bremsstrahlung 
gluon jets originated by the initial quark (the dominant vertex in e + e~ annihilation is 
q — > q + g) and are generated very probably at quite low virtualities (this consideration 
explains the high number of clans and the relatively low population of partons [particles] 
per clan in this case). This interpretation should be compared with what happens very 
reasonably in the same picture in hadron-hadron collisions, where bremsstrahlung gluon 
jets are thought to be generated at quite high virtualities and to have a lot of virtuality 
space for generating larger partonic cascades (the process is dominated here by the gluon 
self- interaction vertex and stronger colour exchanges). 

More accurate analyses of final charged particles MD's at higher energies (at LEP Q, 
UA5 H) revealed violations of the regularity. A shoulder was seen in the MD's both in 
full phase space and in rapidity intervals in e + e~ annihilation and hh collisions, which 
were not described by a single NB(Pascal)MD. The death of the regularity was celebrated 
as an expected and sound fact. Experimental complexity was winning over theoretical 
simplicity. 

A different school of thought pointed out that the NB (Pascal) regularity was not dead 
in multiparticle dynamics, but simply that it was working at a more fundamental level 
of investigation, i.e., at the level of different substructures characterising various classes 
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of collisions 0. The violation of NB (Pascal) regularity in the full sample of events in a 
high energy collision was considered as the indication of the existence of substructures (or 
different classes of events), and the suggestion was to explore the validity of the regularity 
in these substructures. It was shown § that the regularity was surviving in the separate 
2- and 3-jet samples of events in e + e~ annihilation and, presumably, in soft and semi-hard 
components in hh collisions. 

It turns out in fact that the weighted superposition of the mentioned substructures, 
each described by a NB(Pascal) MD, reproduces approximately three observed behaviours 
in final charged particles MD's in both classes of collisions: the first one is the mentioned 
shoulder effect in P n vs n (P n is the MD) [0, ||; the second one are n-oscillations in 
H n = K n /F n vs n fH, |l(J (-^n are factorial moments, K n cumulant moments), and the third 
one is the general behaviour of the forward-backward multiplicity correlation strength 
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The fact that three important effects on which QCD has not yet satisfactory predic- 
tions (the only claim is the onset of the hard gluon vertex) and which can be interpreted 
in terms of the same cause, i.e., the validity of the NB (Pascal) regularity at a more 
elementary level of investigation than initially thought, raises intriguing questions on clan 
properties in all classes of collisions, the main one being: are clans observable, or is clan 
concept a purely statistical one like cluster expansion in statistical mechanics? In or- 
der to approach the problem, we decided to proceed by asking ourselves the following 
preliminary questions. 

• Are clans massive objects? 

• If clans are massive, are clan masses different in different classes of events (or sub- 
structures) in a given collision? 

• If clans are massive, what about clan masses in different classes of collisions? 



2 First question. 

A quantitative answer to the first question has been given by A.Bialas and A. Szczerba fl3|| ; 
it was stimulated by the observed qualitative properties of clan structure analysis when 
applied to multiplicity distributions in rapidity intervals in hadron hadron collisions for 
the full sample of events collected by UA5 Collaboration. Here as discussed previously, the 
charged particle multiplicity distributions in rapidity intervals are of course of NB (Pascal) 
type, P^ B (k,n), with k and n increasing with rapidity and with particles generated by 
each independently produced clan according to a logarithmic distribution, P^(P), with 
(3 = n/(n + k). 

Two assumptions were at the basis of the mentioned generalisation of standard clan 
structure analysis to rapidity intervals: they concern the distributions of clans in rapidity 
and the angular distribution in clan decay respectively. 

As previously mentioned, clans are Poissonianly distributed and independently emitted 
in bremsstrahlung-like fashion. Using energy and (longitudinal) momentum conservation, 
the single-clan (pseudo)-rapidity density has been written |14| as 

^ = A(l-x+) A (l-x_) A (2) 
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with 



m 
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where m is the (average) transverse mass of the clan (uit = \/m 2 + p\, it will be called 
m in following), A is a parameter closely related to the plateau height, that is, to the 
average number of clans per unit (pseudo)-rapidity y, and is the cm. energy. Notice 
that Eq. (H) limits clan emission to the interval \y\ < ln(-y/i/m). 

Assuming furthermore that each clan produces particles according to a logarithmic 
MD, whose generating function is 



9iog{z) 



^ Z nW) ln(l-/3) 
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with the average multiplicity per clan, n c , given by 



{(3 - 1) ln(l - (3) ' 



(5) 



then the generating function for the MD in the (pseudo)-rapidity interval Arj turns out 
to be 111 



G(z; Arj) = exp 



dN 



In 



l-JL. p (y-Arj){z-l) 



dy 



ln(l - 0) 



dy 
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p(y; Arj) is the fraction of particles, produced by a clan of (pseudo)-rapidity y, falling 
within Arj; it was also assumed that for a fixed clan multiplicity the MD of particles falling 
within Ar] is binomial, i.e., particles emitted by each clan are emitted independently from 
each other. The integration is over the full range allowed by the kinematical limits of 
Eq. (§). The Authors of Ref. |13| assumed further, for the probability density function 
within a clan to produce a particle at 7], given that the clan is at y, a form based on the 
hypothesis of isotropic decay: 
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{uj = 1 gives isotropic decay if r\ is pseudo-rapidity; the width of the distribution is 
proportional to uj) and thus computed 



p(y; Ar]) 



Hv,y) dr)- 



At/ 



There are, in summary, 4 free parameters to be used to fit the experimental data: (3, 
A, uj and m. 

Experimental data on pp collisions at 546 GeV are approximately reproduced with the 
following choice of the parameters in the generalised model of Ref. |13| 



A = 0.855, 



m 



3.15, 



UJ 



1.45, = 0.90. 



(9) 



The obtained multiplicity distribution is not indeed of NB type except in full phase 
space (deviations are significant forn < 3, and in k parameter rapidity dependence at the 
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Table 1: Parameters A, m, w and /3, obtained from fitting Eq. (|6|) to data for each component at various 
cm. energies and in different classes of collisions, are shown in the top part; other quantities derived from 

the fit parameters are show in the bottom part. 





(a) 




(b) 


(c) 






pp collisions 


pp 


collisions 


e + e annihilation 




63 GeV 


900 GeV 


91 GeV 




soft 


soft 


semi-hard 


2-jet 


3- jet 


A 


1.14 


0.92 


2.09 


1.60 


3.33 


m (GeV/c 2 ) 


1.80 


1.47 


3.43 


0.62 


1.10 


UJ 


0.84 


1.95 


1.35 


1.34 


0.56 


P 


0.79 


0.83 


0.87 


0.62 


0.59 


dN/dy\ y=0 


1.15 


0.92 


2.06 


1.57 


3.07 


N 


5.59 


9.83 


16.8 


11.0 


16.6 


n c 


2.41 


2.69 


3.21 


1.64 


1.62 



border of phase space). A later analysis of these data || has shown indeed that a shoulder 
effect occurs in the multiplicity distribution and that one should try to explain observed 
data in n-charged particle multiplicity distributions in rapidity in terms of the weighted 
superposition of two classes of events [soft (no mini-jets) semi-hard (with mini-jets)] each 
described by a MD of NB (Pascal) type. This consideration notwithstanding, the study in 
Ref. [ 13j is, in our opinion, very instructive because it allows to determine characteristic 
clan parameters, like its width and mass, under simple assumptions and it can be used to 
determine important properties of the same parameters in different classes of events and 
in different collisions. 

The answer to the first question is therefore positive. 



3 Second question. 

This result led us to the next question. We studied first pp collisions at 63 GeV in rapidity 
intervals: according to our knowledge, at such cm. energy only one component (the soft 
one) is usually assumed to control the dynamics of the collision and the shoulder effect 
is, to a good approximation, negligible. 

The average number of particles per clan, the average number of clans, the dispersion, 
the parameter k and charged particle multiplicities in pseudo-rapidity intervals r] c < 2.5 
are quite well reproduced by the set of parameters shown in Table [TJa: they are obtained 
by fitting, with the least square method, the average multiplicity (n) and the quantity 
D 2 /(n) 2 — 1/ (n) of the distribution obtained from the generating function (^), in terms of 
the four parameters A, m, uj and /3, to the corresponding moments of the NBMD (namely 
n and k" 1 , respectively) for pseudo-rapidity intervals At] = [— r] c , t] c ] with r\ c < 2.5 (see 
Fig-©. 

We remark that 

a. soft events only are considered in pp at 63 GeV cm. energy (in fact, one component 
only determines reasonably well the n-charged particle multiplicity distribution in rapidity 
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Figure 1: Average multiplicity, dispersion, average number of clans and average number of particles per 
clan in different pseudo-rapidity intervals [— y c , y c ] for the soft component of the MD in pp collisions at 
63 GeV: open circles are data and the solid line our fit. 



intervals) ; 

b. the average number of particles per clan, n c , is not correctly reproduced for 
rapidity intervals larger than r] c = 2.5; a clear bending of n c is visible in the data and 
not its increase with the rapidity interval considered as shown in the generalised model 
(clans in central rapidity regions are larger than in more peripheral regions at the border 
of phase space 

This fact has important consequences on the determination of the leakage parameter 
introduced [JXTJ in the study of particles generated from clans lying in one hemisphere 
to the opposite one in forward-backward multiplicity correlations and suggests that the 
leakage parameter should be larger in broader clans than in smaller ones (and not the 
same throughout all the allowed rapidity range as done in [ll|]). The proposed value of the 
leakage parameter in the just mentioned reference should therefore be considered rather 
an average value between leakage parameters of large and small clans than a rapidity 
independent value. 

We decided then to study pp collisions at 900 GeV cm. energy. According to our 
experience here not only one but two components are controlling the dynamics of the pro- 
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Figure 2: Fit to the average multiplicity and dispersion in different pseudo-rapidity intervals [— y c , y c ] 
for the two components of the MD in pp collisions at 900 GeV: soft component: open circles (data) and 
solid line (fit); semi-hard component: filled circles (data) and dashed line (fit). Data points are from 

Ref. @. ' 



cess, i.e., the soft and the semi-hard component [IfJ. We fit therefore the four parameters 
A, u), m and (3, with the method previously explained, for each component separately 
using available UA5 data ||, i.e., the NB fits in pseudo-rapidity intervals [— T] c , t/ c ] with 
i] c = 1 . . . 5; because, as discussed in clans emitted close to their kinematical limit 
ln(-y/i/m) appear to be smaller, we do not make use of full phase space values. The fits 
turn out to be good, as shown in Fig. ^|, except that they do not reproduce well the de- 
crease in the average number of particles per clan close to full phase space. The resulting 
values of the fit parameters are shown in Table [TJ6. Figure |3| shows the densities dN / dy, 
Eq. (§), and 0(7/; 0), Eq. (|7|), for the two components at 900 GeV separately, using the 
results of the fit. 

In comparing parameters behaviour for the soft component at 63 GeV and 900 GeV 
we notice that clan masses and distribution widths vary with cm. energy and the plateau 
height is slowly decreasing. The decrease of the clan mass in the soft component with 
increasing energy could be due to the intentionally overlooked contamination of semi-hard 
events at 63 GeV, as will be discussed in the last section. 
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Figure 3: Clan density dN/dy, Eq. (|^), and single particle pseudo-rapidity probability density in a clan 
0(y;O), Eq. (fj), for the soft (solid line) and semi-hard (dashed line) component at 900 GeV cm. energy; 

parameters from Table [l]&. 

In comparing, next, parameters behaviour for the soft and semi-hard components at 
the same cm. energy (900 GeV) we remark that clan masses and plateau heights are 
much higher in the semi-hard than in the soft component, whereas the distribution width 
is much higher in the soft than in the semi-hard component. This fact shows that heavier 
particles are produced more in semi-hard than in the soft component. 

Interestingly the average number of particles per clan is bending in larger rapidity 
intervals both in the soft and in the semi-hard component suggesting also here that clans 
are larger in central rapidity intervals than in the peripheral ones. 

Accordingly leakage parameters in forward backward multiplicity correlations should 
be larger when clans have larger masses and their particle content is distributed in more 
central rapidity intervals. 

In conclusion clan masses on the average are different in different classes of events or 
substructures of a given class of collision. 

4 Third question. 

It is interesting to check clan properties also in e + e~ annihilation at LEP cm. energy 
within the generalised model of Bialas and Szczerba. 

We apply it to the 2-jet and 3-jet samples of events and we find that data are approx- 
imately reproduced with the choice of parameters shown in Table p]. 

We conclude in comparing the average masses of 2- and 3-jet samples of events in e + e~ 
annihilation at LEP energy with the average mass of the soft component at 63 GeV in 
hh collisions that clan masses in 2-jet and 3-jet events are much lower than in the soft 
component at ISR energies. 

5 Remarks and questions for future experiments 

Assuming the validity of NB (Pascal) regularity and its interpretation in terms of clan 
structure for each substructure (or component or subsample of events) characterising dif- 
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ferent classes of high energy collisions it was shown in previous works that the weighted 
superposition of NB (Pascal) MD's (one for each component) describes quite well exper- 
imental final charged particle MD properties which a single NB (Pascal) MD is unable 
to reproduce. It should be pointed out that QCD has (up to now) quite poor prediction 
on all just mentioned experimental facts: experts claim that they are all consequences 
of the onset of hard gluon radiation and rely fully on purely complicated higher order 
perturbative calculations to be eventually performed in the future. From a stricter, and 
complementary, phenomenological point of view, the above mentioned successes in de- 
scribing data in terms of the weighted superposition mechanism of two (or eventually 
more) NB (Pascal) MD's and the related interpretation in terms of clans of the differ- 
ent substructures of a collision raise intriguing problems, which we summarised in the 
following question: are clans observable objects? 

This paper is an attempt to answer this interesting question. It concerns mainly energy 
and rapidity dependence of average clan masses of the substructures characterising e + e~ 
annihilation and hh collisions. Results are not inconsistent with general expectations and 
support the idea that clans could be indeed observable. 

Some warnings are needed. Our work is based on informations and extrapolations 
not coming from dedicated experiments on the subject, as these are not available at 
present. We are therefore aware of the intrinsic limitations of our approach. For instance 
the separation between soft and semi-hard events at 900 GeV in hh collisions comes (as 
already pointed out) from a fit proposed in order to reproduce observed experimental 
data on MD's. The question still on the carpet is indeed what is a soft and what is a semi- 
hard event in hh collisions. Some progress has been made by CDF Collaboration |I7J , but 
a definite generally accepted answer is lacking. Similarly at 63 GeV cm. energy (ISR) in 
view of the lack of a clean separation between soft and semi-hard events, all events were 
taken to be soft (and consequently a single NB (Pascal) MD was used for describing the 
full sample of events). It cannot be excluded that the relatively high average clan mass 
found in the soft component at 63 GeV cm. energy (m « 1.80 GeV/c 2 ) (no semi-hard 
events were assumed to contribute to the full sample of events) with respect to that of the 
average clan mass of the soft component at 900 GeV (m m 1.47 GeV/c 2 ) (here both soft 
and semi-hard events [according to Fuglesang's fit] were assumed to contribute to the full 
sample of events) might be the consequence of the contamination of a certain percentage 
of semi-hard events at 63 GeV which at present we are unable to disentangle and which 
will modify the calculated average clan mass of the soft component. 

In addition, the separation between 2-jet and 3-jet events in e + e~ annihilation is 
jet-finder algorithm dependent and it might be that one should consider the 3-jet sample 
events separated into two extra subsamples of events (e.g., Mercedes-like events and hard- 
gluon events). All these remarks make of course our conclusions questionable. These 
consideration notwithstanding, it seems to us interesting to explore the possibility that 
clans do have mass. This search should be done in our opinion in future experiments in 
two steps. 

Firstly, it is important to check with dedicated experimental work existing substruc- 
tures in the full sample of events of each class of collisions and verify then that they 
agree with NB (Pascal) behaviour or at least with infinitely divisible distribution (IDD) 
behaviour (which maintains the Poissonian nature of the first step of the production pro- 
cess) . Secondly, one should apply the analysis of the present approach to each component 
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in order to determine clan masses properties. In this respect it is relevant to remember 
that the average masses we are referring to are average transverse masses and therefore 
one should pay attention to the transverse momentum component fraction contributing 
to the real average clan mass which might be disturbing in the separation between soft 
and semi-hard events. 

It is clear that if clans are massive and their average masses vary with energy and 
rapidity the next question is: do clans have other other quantum numbers? If clans 
are real objects, do they continue to be independently produced and to decay with no 
correlations among particles generated by different ancestors at all cm. energies, or there 
exists a threshold beyond which clans loose their original independence, start to overlap 
and then to interact among themselves, modifying the initial simplicity of clan structure 
interpretation of NB (Pascal), or IDD, regularities? 

Other interesting remarks and subtle questions can be added to the previous ones. 
They are too premature in our opinion. The point we want to make, following the results 
of this approach, is that if clans are massive objects (with eventually other quantum num- 
bers) and their masses vary with cm. energy and rapidity, and clan masses are different in 
different substructures of high energy collisions, the hadronization process itself should be 
consistent with the new scenario which emphasises the role of clans with respect to that of 
final charged particles. The hadronization process in this perspective will be dominated 
by massive clans (primaries?) whose search and characteristic properties study will be a 
possible new frontier in multiparticle dynamics in all high energy collisions, a challenging 
problem to experimentalists of the next generation accelerators. 
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